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Abstract 

Tine uncertainty on tiie calorimeter energy response to jets of particles is derived for the ATLAS 
experiment at the Large Hadron Collider (LHC). First, the calorimeter response to single isolated 
charged hadrons is measured and compared to the Monte Carlo simulation using proton-proton col- 
lisions at centre-of-mass energies of v^ = 900 GeV and 7 TeV collected during 2009 and 2010. 
Then, using the decay of K^ and A particles, the calorimeter response to specific types of particles 
(positively and negatively charged pions, protons, and anti-protons) is measured and compared to 
the Monte Carlo predictions. Finally, the jet energy scale uncertainty is determined by propagating 
the response uncertainty for single charged and neutral particles to jets. The response uncertainty is 
2-5% for central isolated hadrons and 1-3% for the final calorimeter jet energy scale. 
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Abstract The uncertainty on the calorimeter energy 
response to jets of particles is derived for the ATLAS 
experiment at the Large Hadron Collider (LHC). First, 
the calorimeter response to single isolated charged hadrons 
is measured and compared to the Monte Carlo simu- 
lation using proton-proton collisions at ccntre-of-mass 
energies of ^/s = 900 GeV and 7 TeV collected dur- 
ing 2009 and 2010. Then, using the decay of Kg and 
A particles, the calorimeter response to specific types 
of particles (positively and negatively charged pions, 
protons, and anti-protons) is measured and compared 
to the Monte Carlo predictions. Finally, the jet energy 
scale uncertainty is determined by propagating the re- 
sponse uncertainty for single charged and neutral parti- 
cles to jets. The response uncertainty is 2-5% for central 
isolated hadrons and 1-3% for the final calorimeter jet 
energy scale. 

Keywords Calorimeter • jet • ATLAS 



1 Introduction 

Partons scattered in proton-proton interactions are mea- 
sured with the ATLAS detector as coUimated jets of 
hadrons. The uncertainty on the jet energy scale is 
the largest source of detector-related systematic un- 
certainty for many physics analyses carried out by the 
ATLAS Collaboration, from the di-jct cross-section and 
top mass measurements, to searches for new physics 
with jets in the final state. It is thus the subject of an 
extensive and detailed study [1]. 

The jet energy measured by the calorimeter is cor- 
rected for calorimeter non-compensation and energy loss 
in dead material. The corresponding jet energy scale 
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correction factor is referred to as the JES. The JES is 
derived from Monte Carlo (MC) simulations by compar- 
ing the calorimeter energy of an isolated reconstructed 
jet to that of the particle jet^ that points to it [1]. 

The uncertainty on the calorimeter energy response 
is a significant component of the total uncertainty on 
the JES. It is derived in this paper by convolving the 
measured uncertainty on the single charged hadron en- 
ergy response and the estimated uncertainty on the neu- 
tral particle energy response with the expected particle 
spectrum within a jet. 

The calorimeter response to single isolated charged 
hadrons, and the accuracy of its Monte Carlo simula- 
tion description, can be evaluated from the ratio of the 
calorimeter energy E to the associated isolated track 
momentum p. The aim of the measurement is to esti- 
mate the systematic uncertainty on jet calorimeter re- 
sponse and therefore the focus is on data-to-MC com- 
parison. The ratio E /p is measured using proton-proton 
collisions at centre-of-mass energies of ^/s = 900 GeV 
and 7 TeV over a wide range of track momenta in the 
central region of the calorimeter. Possible additional 
uncertainties introduced by certain particle species are 
addressed by measuring the response to hadrons iden- 
tified through the reconstruction of known short-lived 
particles. 

This paper is organised as follows. The relevant fea- 
tures of the ATLAS detector are summarised in Sec- 
tion 2. Then the measurement of the single charged 
hadron response, at centre-of-mass energies of -y/s = 
900 GeV and 7 TeV, is discussed in Section 3. Section 4 

^ Particle jets (in MC simulated events) are defined as the 
jets obtained by running the jet finding algorithm on the 
stable particles from the event generator, including those 
with lifetimes longer than 10 ps and excluding neutrinos and 
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describes the measurement of E/p for charged particles 
identified using Kg and A particle decays. In Section 5, 
the determination of the JES uncertainty is discussed. 
This includes the propagation of the charged particle 
response uncertainty to jets and additional systematic 
uncertainties related to the calorimeter response of neu- 
tral particles. The total uncertainty on the calorimeter 
energy response to jets and its correlations are shown 
for jets with transverse momenta between 15 GeV and 
2.5 TeV and in the pseudorapidity^ range \ri\ < 0.8. The 
total JES uncertainty, including all calorimeter regions 
and effects not related to the calorimeter response, is 
discussed in Ref. [1]. Throughout the paper, all uncer- 
tainties are combined in quadrature, unless stated oth- 



covering 1.5 < |?7| < 3.2) are LAr sampling calorime- 
ters with copper absorbers. The copper/tungsten-LAr 
forward calorimeters (FCal) provide both electromag- 
netic and hadronic energy measurements and extend 
the coverage to \r]\ < 4.9. 

The data used for this analysis were triggered us- 
ing the minimum bias trigger scintillators (MBTS) [3]. 
The MBTS are mounted at each end of the detector 
in front of the LAr end-cap calorimeter cryostats at 
z = ±3.56 m and are segmented into eight sectors in 
azimuth and two rings in pseudorapidity (2.09 < \r]\ < 
2.82 and 2.82 < |r;| < 3.84). 



3 Single particle response for charged hadrons 



2 The ATLAS detector 

The ATLAS detector covers almost the whole solid an- 
gle around the collision point with layers of tracking 
detectors, calorimeters and muon chambers and is de- 
scribed in detail in Ref. [2]. Here, the features relevant 
for this analysis are summarised. 

The inner detector (ID) is immersed in a 2 T axial 
magnetic field and provides tracking for charged parti- 
cles with |?7| < 2.5. The ID consists of a silicon pixel 
tracker and silicon microstrip tracker (SCT) covering 
1 77 1 < 2.5 and a transition radiation tracker (TRT) cov- 
ering |?7| < 2.0. 

The calorimeter system covers |?7| < 4.9, using a 
variety of technologies. High granularity liquid-argon 
(LAr) electromagnetic (EM) sampling calorimeters, with 
excellent performance in terms of energy and position 
resolution, cover |?/| < 3.2. They use accordion-shaped 
electrodes and lead absorbers and consist of a barrel 
(EMB, I77I < 1.475) and an end-cap (EMEC, 1.375 < 
|?7| < 3.2). They are longitudinally segmented in depth 
into three layers, with a pre-sampler behind the solenoid. 
For \ri\ < 1.7 hadronic calorimetry is provided by a sam- 
pling calorimeter made of iron and scintillating tiles 
(TileCal). TileCal comprises a large barrel (I77I < 0.8) 
and two smaller extended barrel cylinders (0.8 < \ri\ < 
1.7). It is segmented longitudinally into three layers, 
with a total thickness of about eight interaction lengths 
at ?7 := 0. The hadronic end-cap calorimeters (HEC, 

^ ATLAS uses a right-handed coordinate system with its 
origin at the nominal interaction point (IP) in the centre of 
the detector and the z-axis along the beam pipe. The 3;-axis 
points from the IP to the centre of the LHC ring, and the 
y-eods points upward. Cylindrical coordinates (r, </>) are used 
in the transverse plane, </) being the azimuthal angle around 
the beam pipe. The pseudorapidity is defined in terms of the 
polar angle 9 as 77 = — lntan(S/2). 



In this section, the response of the calorimeters to iso- 
lated charged hadrons is compared to the predictions 
from MC simulation. The ratio of the energy, E, de- 
posited by an isolated charged particle in the calorime- 
ter to the track momentum, p, is used to evaluate the 
uncertainty on the calorimeter response modelling in 
the MC simulation. 



3.1 Event selection 

Events are required to have at least one reconstructed 
vertex with at least four associated tracks. The to- 
tal number of events satisfying this selection is about 
25 million collected in 2009 (approximately one million 
events at y/s = 900 GeV) and 2010 (approximately 24 
million events at y/s = 7 TeV) . 

For every selected event, each track candidate is ex- 
trapolated to the second longitudinal layer of the EM 
calorimeter. A track is defined as isolated if its impact 

point has a distance AR = y {A-qf + {A(t)f > 0.4 
from all other track candidate impact points. The iso- 
lation conditions are studied in detail in Section 4.5. 
The isolated tracks must also have: 

— a transverse momentum of px > 500 MeV"^ , 

— a minimum of one hit in the pixel detector and six 
hits in the SCT, and 

— small transverse and longitudinal impact parame- 
ters (defined in Ref. [4]) computed with respect to 
the primary vertex^, |(io| < 1.5 mm and \zq\ sin 9 < 
1.5 mm. 

^ Below 500 MeV a charged particle loops in the ID and 
does not reach the barrel calorimeter. 

■* In case of multiple vertices, the primary vertex is taken 
to be the one for which the sum of the square of momenta of 
the attached tracks p^ is the largest. 
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The above requirements ensure a good quality of 
the track and reduce contributions from fake tracks to 
a neghgible level. 



reconstructed track momentum. Note that because of 
calorimeter noise fluctuations, E (and therefore E/p) 
can assume negative values. 



3.2 Monte Carlo simulation 

A sample of about 10 (20) million non-diffractive proton- 
proton collision events at ^/s = 900 GeV {^/s = 7 TeV) 
are generated using Pythia 6.421 [5] with the ATLAS 
minimum bias tune 1 (AMBTl) [4]. For isolated, high- 
momentum (pt > 15 GeV) tracks this corresponds to 
^ 60% of the available events in data. All the events 
are run through a full detector simulation [6] based on 
Geant4 [7] . The set of Geant4 physics models used is 
QGSPJ3ERT [8]. The reconstruction and analysis soft- 
ware used for the MC simulation is the same as for the 
data. 



3.3 Defining the E/p observable 

The sum of the energy deposits in layers of calorime- 
ter cells associated to a selected track is computed us- 
ing topological clusters [9] at the electromagnetic scale, 
i.e., without applying any correction for the calorimeter 
non-compensation or for energy loss in dead material. 
The purpose of the topological clustering algorithm is 
to identify areas of connected energy deposits in the 
calorimeter, based on the significance of the energy de- 
posits in cells with respect to the expected noise level. 
Topological clusters are formed around cells with en- 
ergy |£'ccii| > 4o-„oisc ("seeds"), where tTnoiso is the RMS 
of that cell noise. Then, iteratively, the cluster is ex- 
panded by adding all neighbouring cells with li^ceul > 
2o'noisc- Finally, the cells surrounding the resulting clus- 
ter are added, regardless of their energy. The r] — cj) 
position of a cluster z in a given calorimeter layer j, 
i%i ' ^c1 ) ^^ computed as the energy-weighted position 
of the cells in layer j belonging to the cluster. 

The position of the track k extrapolated to the layer 
j is {%^ , 0t j? ) . The energy of a cluster in the layer j [Ej ) 
is associated to the track if: 



(%r -r]ciy + 



b'l^ 



K'l)' 



<R 



coll- 



(1) 



The parameter i?coii is set to 0.2 based on a trade- 
off between maximising the particle shower contain- 
ment and minimising the background contribution com- 
ing from neutral particles produced close to the track. 
Roughly 90% of the shower energy is collected in a cone 
of such size. 

The energy E associated to a track is computed as 
the sum of the associated energy deposits in all lay- 
ers, E ^ ^ Ej, and the ratio E/p is formed with the 



3.4 E/p distributions 

The E/p distributions in two representative regions of 
rj and track momentum are shown in Fig. 1. The large 
number of entries with E/p = corresponds to isolated 
tracks that have no associated cluster in the calorime- 
ter. Several effects may be responsible for this: 

— Particles can interact hadronically before reaching 
the calorimeter (in the ID, cryostat or solenoid mag- 
net). Such particles can change their direction, or 
produce a large number of low momentum secondary 
particles. 

— A cluster is created only if a seed is found. Hadrons 
with low momentum and an extended shower topol- 
ogy sometimes do not have a single cell energy de- 
posit large enough to seed a topological cluster. 

The cases where the calorimeter response is com- 
patible with zero have been further studied. The prob- 
ability that the calorimeter response is suppressed by 
noise threshold requirements, P{E = 0), is shown in 
Fig. 2a as a function of the amount of material (in nu- 
clear interaction lengths) in front of the active volume 
in the central (|?7| < 1.0) calorimeter region. When the 
hadron passes through more material, the probability 
that no energy is associated to the track increases. 

Figure 2b shows P{E = 0) as a function of the 
track momentum in the central region of the calorime- 
ter {\r]\ < 0.6). In this region, the dead material in front 
of the calorimeter is approximately constant. The prob- 
ability decreases with increasing track momentum. In 
general, P{E = 0) is well predicted by the MC simula- 
tion. 



3.5 Background subtraction 

The energy measured inside the cone of AR < i?coii = 
0.2 centered around the track impact point may be 
contaminated by energy deposits from the showers of 
close-by particles produced in the proton-proton colli- 
sion. The track isolation requirement suppresses possi- 
ble shower contamination from charged particles. There 
is no obvious way to suppress shower contamination 
from photons, mostly produced in tt'^ — )■ 77 decays, 
and neutral hadrons. The neutral particle background 
contribution to the E/p measurement in the MC sim- 
ulation depends on the event generator settings of the 
parameters governing non-pcrturbativc QCD processes 
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Fig. 1: (a) E/p distribution for isolated tracks "with an impact point in the region \t]\ < 0.6 and "with a momentum 
in the range 1.2 < p < 1.8 GeV. (b) E/p distribution for tracks "with impact points in the region 1.9 < |?7| < 2.3 
and "with momenta in the range 2.8 < p < 3.6 GeV. 
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Fig. 2: (a) Probabihty to measure a caforimeter response consistent -with zero P{E = 0) as a function of the 
amount of material in nuclear interaction lengths in front of active volume of the calorimeter for \r]\ < 1.0. (b) 
P{E = 0) as a function of track momentum for I77I < 0.6. 



and on the modelling of the calorimeter response to lo"w 
momentum neutral particles, and it is therefore diffi- 
cult to model correctly. The neutral background is thus 
subtracted from the measured response using an in situ 
background estimate. In the follo"wing, the expression 
"EM (HAD) energy" -will refer to the energy deposited 
in the electromagnetic (hadronic) calorimeter. 

The background subtraction relies on the assump- 
tion that the EM energy from photons and neutral 
hadrons is independent of the energy deposited by the 
selected track. Charged hadrons are selected that be- 
have like minimum ionising particles in the EM calorime- 
ter and start their sho"wer in the hadronic calorimeter 



(late-sho"wering hadrons). Excluding a narro"w region 
around the late-sho"wering hadron track, the remain- 
ing EM energy is mainly due to sho"wers from neutral 
particles. The strategy is sketched in Fig. 3. 



Late-sho"wering hadrons are selected by requiring a 
small amount of EM energy in a cone of AR < 0.1, 
-^EM < 1-1 GeV, and a large HAD energy fraction, 
-^had/p > 0-4- The background is measured in the EM 
calorimeter in an annulus around the late-sho"wering 
charged hadrons. The mean of the background distribu- 
tion over many events in a given momentum and pseu- 
dorapidity bin estimates the energy deposition of pho- 
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'^(ICH) 



been validated by comparing its prediction with that of 
a more complex procedure, which is described in Ap- 
pendix A. 

The background subtraction procedure is applied in 
the following to all E/p measurements, both at ^/s = 
900 GeV and ^/s = 7 TeV. Because the background con- 
tribution increases with ^/s, the JES uncertainty esti- 
mation relies on the 900 GeV data for the measurement 
of calorimeter response to tracks with momenta below 
2.2 GeV. 



Fig. 3: Sketch of the background subtraction method, 
selecting isolated charged hadrons (ICH) that pass the 
EM calorimeter as minimum ionising particles (MIP) 
and shower in the hadronic calorimeter. 



tons and neutral hadrons showering in the EM calorime- 
ter: 



{E/p) 



BG 



E;0.2 



PO.I 
-^EM 



(2) 



where E^^ is the EM energy in a cone of AR < 0.2. 

The background contribution from neutral hadrons 
depositing their energy in the hadronic calorimeter was 
estimated with a similar technique applied to different 
hadronic calorimeter layers and found to be negligibly 
small. 

Figures 4 and 5 show {E/p)-qq as a function of p in 
two bins of pseudorapidity at ^/s = 900 GeV and ^/s = 
7 TeV. The background {E/p)bg at ^ = 7 TcV in 
both bins is about 0.04-0.08 for p<W GeV (0.02-0.08 
at ^/s = 900 GeV) and decreases to ^ 0.04 for track 
momenta of ^20 GeV. The general trend is confirmed 
by the MC simulation, although the two show some 
significant differences. The discrepancy is attributed to 
an imperfect modelling of non-pcrturbative QCD pro- 
cesses by the Pythia settings used for the MC event 
simulation. 

The quantity that is studied in detail in the next sec- 
tion as a function of the track momentum and calorime- 
ter impact point pseudorapidity is: 



(E/p) = (E/p),,^ - ^{E/p)bg 



(3) 



where (i*^/p)raw is the mean value of the E/p distribu- 
tion before background subtraction (cf. Fig. 1). 

The background (Eq. 2) is rescaled by the factor of 
4/3, which is the ratio of the area of the full AR < 
0.2 cone to that of an annulus with 0.1 < AR < 0.2, 
to take into account the background contribution in 
^EM J since the background is only measured in this an- 
nulus. A uniform energy density of the background in 
the AR < 0.2 cone is assumed. The assumption has 



3.6 Results 

Several systematic uncertainties on the measurement 
are estimated. Each is taken to be completely correlated 
between all pseudorapidity and momentum bins in the 
{E/p) measurement: 

— Track selection: The dependence of {E/p) on the 
track selection requirements in Section 3.1 has been 
estimated by varying the number of silicon tracker 
hits required and the impact parameter selection 
with respect to the primary vertex within reasonable 
ranges. The MC-to-data ratio of {E/p) was found to 
be almost unaffected by variations in the track se- 
lection. The maximum variation found in the ratio 
(0.5%) is taken as a systematic uncertainty. 

— Track momentum scale: The uncertainty on the mo- 
mentum scale p as measured by the inner detector is 
negligibly small for p < 5 GeV [10]. For p > 5 GeV, 
a conservative 1% uncertainty has been assumed on 
the momentum scale. 

— Background subtraction: The difference between the 
background estimate obtained with the method de- 
scribed in Section 3.5 and with the validation method 
described in Appendix A is taken as a systematic 
uncertainty. This results in a 1% uncertainty on the 
E/p measurements. 

The mean E/p value after background subtraction 
is evaluated in bins of momentum and pseudorapidity. 
Figures 6 and 7 show {E/p) as a function of the track 
momentum, in two different \'q\ bins up to \'q\ ~ 1.1, at 
yi = 900 GeV and y^ = 7 TeV. The lower parts of the 
figures present the ratio of MC simulation to data. The 
data with ^/s = 900 GeV and ^/s = 7 TeV agree within 
the statistical uncertainty. The maximum momentum 
that can be probed with the data considered is approx- 
imately 30 GeV. The agreement between data and MG 
simulation is within ~2% for particles with momenta in 
the 1-10 GeV range, and it is around 5% for momenta 
in the 10-30 GeV range. Below 1 GeV, where tracks are 
just at the kinematic threshold of entering the calorime- 
ter volume, large differences of ~10% or more between 
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Fig. 4: {E/p)bg as a function of the track momentum at -y/s = 900 GeV for (a) \ri\ < 0.6 and (b) 0.6 < \ri\ < 1.1. 
The black niarkers represent the background estimated from cohision data, while the green rectangles represent 
the MC prediction, with the vertical width showing its statistical uncertainty. The lower panes show the ratio of 
the MC prediction to collision data. The MC prediction of the background energy deposit is obtained using the 
same procedure as applied to the data. 
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Fig. 5: {E/p)bg as a function of the track momentum at y/s — 7 TeV for (a) jryj < 0.6 and (b) 0.6 < \r]\ < 1.1. 
The black markers represent the background estimated from collision data, while the green rectangles represent 
the MC prediction, with the vertical width showing its statistical uncertainty. The lower panes show the ratio of 
the MC prediction to collision data. The MC prediction of the background energy deposit is obtained using the 
same procedure as applied to the data. 



Single hadron response measurement and calorimeter jet energy scale uncertainty 



data and MC simulation are visible. However, due to 
the low absolute calorimeter response to very low mo- 
mentum particles, these differences are not critical for 
the JES determination. 



4 Calorimeter response to identified hadrons 

The extrapolation of the previous single particle re- 
sponse studies into the environment of a jet requires 
understanding of two additional effects. A jet includes 
a variety of hadrons that may differ from the inclu- 
sive sample of isolated hadrons. Therefore, measuring 
the average response to different species of particles 
is valuable to ensure that the Monte Carlo correctly 
models all aspects of the jet shower. Additionally, the 
hadrons in a jet are not isolated. Threshold effects and 
hadronic shower widths affect the calorimeter response 
to the multi-hadron system. In order to address these 
two points, this study adds two additional features to 
help complete the understanding of calorimeter response. 
To minimise the impact of the background and to sim- 
plify its estimation, the measurements are presented as 
a difference between the ratio E/p for two different par- 
ticle types or as a ratio to the inclusive measurement. 

First, single hadrons arc identified using decays of 
Ks (for positive and negative pions), A (for protons), 
and A (for anti-protons) particles. These single hadrons 
are required to be isolated from all other charged par- 
ticles in the event. Single pions will have manifestly 
lower energy response distributions from those of anti- 
protons, because of the eventual annihilation of the 
anti-proton. By identifying and isolating single pions, 
single protons, and single anti-protons, the effects of 
hadronic interactions and annihilation can be separated 
at low to moderate energies, where they are most im- 
portant. 

Second, when the mother particle is highly boosted, 
the decay products are more coUimated. For track mo- 
menta of ~2-6 GeV from Ks decays, the range of open- 
ing angles of the decay products allows a measurement 
of calorimeter response as more tightly coUimated pairs 
of pions are selected. There are two related effects probed 
by such a measurement. Energy deposited in the calorime- 
ter may fall below the thresholds for reconstruction and 
be neglected as consistent with noise. As two showers 
overlap, the addition of energies that might have in- 
dividually been below threshold can produce a signal 
above these noise thresholds. The increase of the signal 
is related to both the thresholds themselves and the 
width of the hadronic shower. For this study, the tt+tt" 
system is required to be isolated from other charged 
particles in the calorimeter. 



4.1 Event selection and observable definition 

The event selection for this measurement follows closely 
the event selection of the inclusive measurement already 
presented, except that a secondary vertex is required in 
the event. The events are collected using random trig- 
gers at y/s = 7 TeV and correspond to an integrated lu- 
minosity of about 800 /ib^^. The same primary vertex 
requirements are applied, and the same requirements 
are placed on tracks entering the measurement of E/p, 
except the requirements on the tracks' impact parame- 
ters. The Monte Carlo used is described in Section 3.2. 
The definition of the ratio E/p and the isolation re- 
quirement for tracks is described in Section 3.3. 

The ratio E/p is presented as a function of avail- 
able energy, Ea- For pions, this is simply the particle's 
total energy: Ea 

m? — m. For 



yjp'^ +m^. For protons, only the 
kinetic energy is included: Ea 



vV 



anti-protons, the kinetic energy plus double the rest- 
mass is included, in order to take into account anni- 
T? + m. The available energy 



hilation, Ea 



V^ 



is therefore calculated using the information from the 
tracker. 



4.2 Reconstruction of short-lived particle candidates 

The reconstruction and selection of long-lived particles 
is based on previous ATLAS results [10]. The decay 
Ks — > tt+tt", which dominates the Ks decays to other 
charged particles by several orders of magnitude, is used 
to identify pions. The decay A — )■ Tr~p {A — > ir'^p), 
which dominates the A (A) decays to other charged par- 
ticles by several orders of magnitude, is used to identify 
protons (anti-protons). Both decay product tracks are 
required to have px > 100 MeV, and the tracks used 
for the E/p measurement must have pt > 500 MeV. In 
the case of A (A) candidate decays, the charge of the 
higher momentum track is used to label the candidate 
as a yl (positively-charged higher momentum track) or 
A (negatively-charged higher momentum track), as is 
kinematically favored. The tracks entering the E/p dis- 
tributions are additionally required to have a pseudorapidity- 
dependent number of hits in the TRT, the outermost 
tracking system. This requirement suppresses the ef- 
fect of nuclear interactions in the material of the inner 
detector, particularly in the outer layers of the silicon 
tracker. 

The individual tracks entering the response distri- 
butions arc divided into two bins of pscudorapidity, 
jryl < 0.6 and 0.6 < \ri\ < 1.1. Several higher pscu- 
dorapidity bins show consistent results, albeit with sig- 
nificantly lower statistics. Example distributions of re- 
constructed mass for candidates with at least one cen- 
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Fig. 6: (E/p) at y/s = 900 GcV as a function of the track momentum for (a) \rj\ < 0.6 and (b) 0.6 < \rj\ < 1.1. 
The black markers represent the colhsion data, while the green rectangles represent the MC prediction, with the 
vertical width showing its statistical uncertainty. The lower panes show the ratio of the MC simulation prediction 
to collision data. The grey band indicates the size of the systematic uncertainty on the measurement. The dotted 
lines are placed at ±5% of unity and at unity. 
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identical central values for both Gaussians, and for the 
A and A, a modified Gaussian is used of the form 
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Fig. 8: (a) The reconstructed mass peak of vl candidates 
with central {\ri\ < 0.6) proton-candidate tracks in data 
(black points) and MC simulation (green histogram), 
(b) The reconstructed mass peak of A candidates in 
data. The distribution is shown separately for positive 
(i.e. proton candiate) and negative (i.e. pion candidate) 
tracks in the central region {\r]\ < 0.6). 



y = a X exp 



-L IJ i 

-- X X 1+-/2 

2 



(4) 



where x = — — , 771 is the reconstructed mass, and a, b 
and c are free parameters of the fit. The peak mean, b, 
is stable to within a few hundred keV over all pseudo- 
rapidity bins. The fitted width of the mass peak, c, in- 
creases with track pseudorapidity, due to the track reso- 
lution. Example fit results for the A are shown in Fig. 8 
for positively charged tracks (proton candidates) and 
negatively charged tracks (pion candidates) in the cen- 
tral pseudorapidity bin for the E/p measurement. The 
negatively charged tracks are softer, which is reflected 
in the smaller statistics for the negatively charged track 
mass distribution. 

For each mass peak and each bin of track pseudora- 
pidity, an acceptance window is constructed in order to 
optimise both signal purity and statistics. The width of 
the window is set at three times the width of the nar- 
rower Gaussian for the Ks and three times the width of 
the modified Gaussian for the A {A). The purities are 
over 97% for the Ks candidates and above 92% for the 
central A candidates. The purities in the MC simula- 
tion are within 2% of those in the data. The number of 
tracks from signal candidates, extracted from the fits, 
in the data and 20 million MC simulation events are 
shown in Table 1. Roughly twice as many Ks candi- 
dates and three to four times more A and A candi- 
dates are found in the data than in the MC simula- 
tion, because of the limited MC simulation statistics. 
The difference in Ks and A yields between data and 
MC simulation are discussed further in Ref. [10], and 
the response measurements here are insensitive to these 
differences. Only the positively and negatively charged 
tracks that form short-lived particle candidates within 
the defined acceptance window are considered in the 
remainder of this section. 



tral (|?7| < 0.6) track are shown in Fig. 8. The mass 
peaks stand out clearly over the background. The com- 
position of the background will be discussed further in 
Section 4.3. 

The signal purity for a given pseudorapidity inter- 
val, where the binning is always in terms of the kine- 
matic properties of the track of interest, is estimated 
using a fit to the signal over a cubic polynomial fit to 
the background, following Ref. [10]. For the Ks, a dou- 
ble Gaussian signal function is used, constrained to have 



Particle 



\rj\ < 0.6 
Data MC 



0.6 < |r;| < 1.1 
Data MC 



7r+ from Ks 
iT~ from Ks 
p from A 
p from A 



4.04x105 2.04x10^ 

3.94x10^ l.gSxlO'^ 

2.63x104 8.03x103 

2.64x10'* 7.21x10^ 



2.52x10^ 1.28x10^ 

2.49x10^ 1.26x10^ 

1.81x10'* 5.62x10^ 

1.53x10* 4.10x103 



Table 1: The number of signal candidate tracks, ex- 
tracted from the fits described in the text, in data and 
MC simulation. 
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4.3 Background subtraction 

There are three sources of charged backgrounds enter- 
ing the E/p distributions. Nuclear interactions in the 
material of the inner detector that fake short-lived par- 
ticle candidates are suppressed by the narrow mass ac- 
ceptance window. The smoothness of the distribution 
of secondary vertex positions confirms that these are at 
the level of a few percent. The tracks may also undergo 
nuclear interactions prior to entering the calorimeter. 
These interactions are suppressed by requiring that the 
track have hits in the TRT. An additional charged back- 
ground comes from combinatorics, particularly in high 
track-multiplicity events. Because the purity of Ks can- 
didates is typically 5-10% higher than that for A ox A 
candidates, and because the majority of the background 
is from charged pions, charged background corrections 
are applied only for calculating proton and anti-proton 
response. The response to pions is taken from the Ks 
candidates and is used to correct the response to pro- 
tons and anti-protons. 

When charged pions enter the proton energy re- 
sponse distribution, they are given an incorrect track 
mass hypothesis. The track extrapolation to the calorime- 
ter does not calculate energy loss differently for pions 
and (anti-)protons, so that the measured associated en- 
ergy in the calorimeter remains the same. However, the 
available energy changes significantly when the differ- 
ent mass hypothesis is used. Thus, in order to subtract 
the charged pion background from the (anti-)proton 
distribution, the pions from Ks candidates are given 
the (anti-)proton mass hypothesis and their response is 
re-calculated. The (anti-)proton response in a bin i of 
available energy and 77, {E /p)i, is then given by 



{E/p), = ^{{E/p)T 



(1 - 60 {e/pY,) 



(5) 



where e^ is the purity of the sample in that pscudora- 
pidity bin, {E/pY{™ is the measured response to the 
(anti-)protons, and (E lp)1 is the measured response to 
the charged pions with a proton mass hypothesis. 

The corrections are small in the region \r]\ < 1.1, 
falling from ^5% at low available energy to ~1% at 
high available energy for both protons and anti-protons. 
The corrections are derived independently for data and 
MC simulation to ensure that the differences in purity 
and any differences in response are taken into account. 

As the momenta of the daughter tracks increase, 
the position resolution of the secondary vertex broad- 
ens. This introduces a larger combinatorial background 
in the high-momentum bins. The purities in MC sim- 
ulation follow, to ^5%, those of the data with increas- 
ing momentum. The systematic uncertainty on the cor- 
rected value of E/p introduced by the pT-depcndence 



of the purity is determined to be well below 1% and is 
therefore neglected. 

There is an additional contribution to the (anti-) 
proton response measurement for A (A) decays faking 
A (A) decays. Particularly when the yl or yl has low 
energy, the pion may be the higher-momentum track. 
These "fakes" are suppressed by the kinematic cuts ap- 
plied in the candidate selection, but are still present at 
some level, particularly in the lowest bin of available en- 
ergy. This background is dominated by well-understood 
two-body decay kinematics. It should be well-described 
by the MC simulation and is therefore taken into ac- 
count using MC predictions. 

As discussed in Section 3.5, there is an additional 
contribution to E/p from neutral particles in the event, 
since isolation is only defined relative to charged parti- 
cles. Only differences in response between particle species 
are reported here, because the neutral background should 
cancel in the difference (except insofar as the Ks or A 
production occurs near additional activity). The charged 
particle isolation criterion should be sufficient to ensure 
that the neutral background is uncorrelated with any 
jet-like activity in the event. The cancellation of the 
neutral background is tested using a simulation sample 
of single particles and is found to be valid up to the 
available statistical accuracy. Therefore, no additional 
correction or uncertainty is added for this background. 
Thus, background systematic uncertainties are found to 
be negligible with respect to the statistical errors of the 
MC simulation sample. 

4.4 Isolated identified single particle response 

The uncorrected distribution of E/p for 7r~ , p, tt"*" and 
p in a single bin of available energy and pseudorapidity, 
2.2 < Ea < 2.8 GeV and |?7| < 0.6, is shown in Fig. 9 (cf. 
Fig. 1 for inclusive hadrons in a different p range). All 
the distributions have a small negative tail from noise in 
the calorimeter and a long positive tail from the neutral 
background. There is a much more prominent positive 
tail in the p response distribution, where annihilation 
plays a significant role. A significant fraction of tracks 
have E = 0, for which no cluster of energy was found 
near the track in the calorimeter, as discussed in Sec- 
tion 3.4. The average response, (E/p), is defined as the 
arithmetic mean of the full distribution. The response 
of protons and positively charged pions is well modeled 
in the MC simulation. The response to pions is signifi- 
cantly lower than that to anti-protons in this available 
energy range, although fewer pions have E — 0. The 
agreement in the fraction of pions with E — builds 
confidence in the modelling of the material in front of 
the calorimeter. 
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In order to reduce the effects of neutral background 
subtraction and focus on the differences in response 
amongst particle species, the difference in response be- 
tween several pairs of particles is measured. Figure 10 
shows the difference in response between tt"*" and 7r~ 
in the central (|?y| < 0.6) and forward (0.6 < \ri\ < 
1.1) pseudorapidity bins. The response to 7r+ is higher 
than that to 7r~ at low available energy, in agreement 
with Ref. [11], wherein this difference is attributed to 
a material dependent charge-exchange effect. The dif- 
ference may also be related to the "Barkas correction" 
described in Ref. [12]. The MC simulation is, however, 
reasonably consistent with the data. 

Figure 11 shows the difference in response between 
7r+ and p in the central {\r]\ < 0.6) and forward (0.6 < 
\ri\ < 1.1) pseudorapidity bins. Again, the difference 
between 7r+ and proton response in MC simulation is 
consistent with the data over the entire range of avail- 
able energies. At low available energy, a larger average 
fraction of the initial hadron energy is converted into an 
electromagnetic shower for pions than for protons [12]. 
This leads to a lower response for protons, particularly 
at low available energy, due to the non-compensation 
of the calorimeter. 

Figure 12 shows the difference in response between 
7r~ and p in the central {\r]\ < 0.6) and forward (0.6 < 
1 77 1 < 1.1) pseudorapidity bins. The difference shows a 
~30%x {E/p) disagreement between data and MC sim- 
ulation for E'a < 3 GcV, though they are consistent 
for Ea > 4 GcV. At these low available energies, the 
anti-proton response should be dominated by the an- 
nihilation and the subsequent shower. The difference 



indicates large contributions from processes not well- 
modelled by Geant4. To test the contribution from 
the calorimeter acceptance to this effect, the response is 
constructed excluding the tracks with E = 0. The pions 
and protons show the same level of agreement between 
data and MC simulation, and the disagreement in the 
response difference between pions and anti-protons re- 
mains near 10%. 



4.5 Response to nearby particles 

To study the effects of calorimeter thresholds and noise 
suppression as more tightly collimated pairs of parti- 
cles are selected, isolation is not required with respect 
to the other daughter of the particle candidate. That 
is, when constructing the energy response for the tt"*" 
from a Ks decay, the tt"*" is required to be isolated from 
all tracks in the event other than the tt^ from the Ks 
decay. The response is then examined as a function of 
the distance between the two pions in rj — (j) after ex- 
trapolation to the second layer of the electromagnetic 
calorimeter. No isolation criteria are applied to the 7r~. 
Because the pions are required to come from a particle 
with a narrow width, the kinematics of the second pion 
are constrained, thus removing one possible source of 
discrepancy between data and MC simulation. 

The ratio E/p as calculated here for a single pion 
should increase as another particle approaches it and 
more of the second particle's energy is included in the 
calorimeter energy, E. As shown in Fig. 13 and 14, the 
response docs not vary with separation for large ex- 
trapolated distances, where the single pion is isolated. 
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Fig. 11: The difference in {E/p) between 7r+ from Ks candidates and p from A candidates in traclcs with (a) 
I77I < 0.6 and (b) 0.6 < \r]\ < 1.1. The responses to the protons are corrected for charged particle backgrounds (see 
text). 



The distributions are normahsed to the average E/p 
in the same r/ and p bin in order to remove any dif- 
ferences from the average response. In Fig. 14, pions 
of both charges arc considered in order to increase the 
available statistics. The response rises below AR w 0.3, 
which confirms that the isolation criterion used in Sec- 
tions 3.1 and 4.4 does not lead to a bias. It increases 
significantly for AR < 0.2, when the second pion is 
inside the cone in which the energy is counted when 
constructing E/p. The height of the peak at low AR is 
related to the kinematics of the Ks decay. In particu- 



lar, some of the low AR and low available energy bins 
contain mostly asymmetric decays. 

The height of the response distribution for close- 
by pions and small opening angles is connected to the 
threshold effects, since it is a measurement of the out-of- 
clustcr energy, to first order. The slope of the response 
curve is closely related to the single-hadron shower width. 
The agreement between data and MC is good over most 
of the range of opening angles, but the peak at low 
opening angles is wider in data, indicating a some- 
what broader shower. This discrepancy is qualitatively 
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Fig. 12: The difference in (E/p) between tt^ from Ks candidates and p from A candidates in tracks with (a) 
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Fig. 13: (E/p) for pions as a function of the extrapolated distance between that pion and a pion of the opposite sign. 
In all cases the pions are daughters of a reconstructed Ks candidate and are required to be isolated from all tracks 
in the event except the other daughter of the decay. The response is shown for low energy (2.2 < Ea < 2.8 GeV) 
(a) TT"*" and (b) 7r~ in the central pscudorapidity bin (|?/| < 0.6). The MC simulation has no tracks passing the 
selection in the smallest bin (smallest two bins) of opening distance for 7r+ (tt"). 



in agreement with the discrepancies observed in lateral 
hadronic shower shape, seen in test-beam studies and 
elsewhere [8,13,14,15]. 



5 Calorimeter jet energy scale uncertainty 

The jet energy scale calibration ( JES) corrects the mea- 
sured jet energy for several effects, including calorime- 



ter non-compensation and energy loss in dead material. 
The calibration itself is derived from MC simulation. 
The calorimeter uncertainty on the JES is calculated 
from the uncertainty on the energy response of all par- 
ticles contributing to a jet. Within the MC simulation, 
the energy contribution of each individual particle to a 
given jet can be separated. The convolution of the un- 
certainty on the single particle energy response with the 
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Fig. 14: (E/p) for pions as a function of the extrapolated distance between that pion and a pion of the opposite 
sign. In all cases the pions are daughters of a reconstructed Ks candidate and are required to be isolated from 
all tracks in the event except the other daughter of the decay. The response is shown for (a) more forward 
(0.6 < I77I < 1.1) pions with 2.2 < Ea < 2.8 GcV and for central (|?7| < 0.6) pions in three bins of available energy: 
(b) 2.8 <Ea< 3.6 GeV, (c) 3.6 < Ea < 4.6 GcV and (d) 4.6 < Ea < 6.0 GcV. The MC simulation has no tracks 
passing the selection in the smallest opening distance bin. 



MC jet particle composition is then used to calculate 
the calorimeter uncertainty on the jet energy scale. 

The calorimeter JES uncertainty is derived for the 
well-understood central region of the calorimeter. The 
main reasons for the restriction to the central calorime- 
ter region are the smaller amount of material in front 
of the calorimeter and the existence of combined test 
beam measurements with a setup very similar to the 
final ATLAS configuration. Therefore, the calorimeter 
JES uncertainty is only evaluated using the single par- 
ticle response for \ri\ < 0.8. The total JES uncertainty, 
including all calorimeter regions and all effects not re- 



lated to the calorimeter energy response, is discussed in 
Ref. [1]. 

The JES uncertainty is determined for jets recon- 
structed from topological clusters [9] with the anti-fct jet 
algorithm [16] implemented in the FASTJET package 
[17] for R = 0.4 and R = 0.6 jet sizes. 

The analysis is performed with inclusive di-jct Monte 
Carlo events simulated with Pythia. Jets are selected 
requiring a separation of AR > 2.0 to any other jet with 
EM scale^ pt(EM) > 7 GeV. The jet transverse mo- 

^ The jet energy at the EM scale is the measured energy 
of a jet before applying any corrections for upstream energy 
losses and the non-compensating nature of the calorimeter. 
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mentum is calibrated from the EM scale to the hadronic 
scale using a MC-based JES [1]. 

The numerical evaluation of the uncertainty on the 
jet energy scale is performed with Monte Carlo pseudo- 
experiments. In each pseudo-experiment, the jet en- 
ergy scale is calculated after randomly changing the 
Monte Carlo single particle energy response within the 
appropriate uncertainty range given by the measured 
data/MC ratio. The final uncertainty on the jet energy 
scale is then given by the spread of the distribution of 
the jet energy scale over all pseudo-experiments. Within 
each pseudo-experiment, all randomly changed param- 
eters are kept fixed, such that the energy response cor- 
relations are properly taken into account. 

The uncertainties on the particle energy response 
functions entering the calculation are taken from the 
E/p measurements described in Sections 3 and 4, ATLAS 



ancies of ~ 5% at low p that disappear within the 
statistical uncertainties for p > 10 GcV. 
Out-of-cone AR > 0.2 energy deposits: statistical 
consistent results were found for the E/p data to 
Monte Carlo ratio for calorimeter energy measured 
in AR < 0.3. Hence no additional uncertainty is 
assumed. 



5.2 Additional uncertainty in the jet energy scale 

The E/p measurements only cover the response of charged 
hadronic particles with momenta less than ^ 20 GeV. 
However, depending on the jet momentum, on average 
between 35% and 90% of the energy in jets is carried by 
particles that are not measured in situ using the isolated 
track analysis (mostly photons from 7r° decays, neutral 



combined test beam (CTB) measurements [13] and GEANT4hadrons and high momentum charged hadrons). Hence, 
Monte Carlo predictions. The details arc described in the uncertainty on the energy response to these parti- 
the following subsections. cles is needed in order to obtain the total calorimeter 

uncertainty on the jet energy scale. 



5.1 Additional E/p uncertainties contributing to the 
jet energy scale 

When the single charged hadron response from E/p is 
used to assess the uncertainty on the jet energy scale, 
further systematic uncertainties that might affect the 
propagation of the response to the jet have to be taken 
into account: 

— E/p acceptance: it was shown in Section 3.3 that the 
probability to find E/p = is strongly correlated to 
the amount of upstream material and hence a good 
measure of the E/p acceptance. A fully correlated 
(in p and rj) 28% uncertainty is derived from the 
maximal observed difference between data and MC 
simulation in this probability. 

— Topological clustering effect: given an amount of en- 
ergy released in the calorimeter, the energy collected 
in the topological clusters may differ from the en- 
ergy released in the calorimeter, depending on how 
isolated the energy deposit is. This can introduce a 
bias in the particle response in a jet with respect to 
the response measured for isolated hadrons. A con- 
servative systematic uncertainty has been computed 
by comparing the results of the {E/p) measurement 
obtained using clusters as described in Section 3 to 
those obtained by repeating the measurement us- 
ing all calorimeter cells in a cone of size AR < 0.2. 
The double ratio of data/MC cluster response to 
data/MC cell response is used to estimate the uncer- 
tainty. The relevant result for the central calorime- 
ter region is presented in Fig. 15, showing discrep- 



5.2.1 High momentum charged particles 

In 2004, an ATLAS Combined Test Beam (CTB) pro- 
gram was carried out at CERN. A "slice" of the ATLAS 
detector composed of the final versions of all sub-detectors 
in the barrel region was exposed to Super Proton Syn- 
chrotron (SPS) test beams. The layout of the sub-detectors 
was designed to be as close to that of ATLAS as pos- 
sible. The setup was used to measure the combined 
calorimeter response to single charged pions of ener- 
gies between 20 and 350 GcV for pseudorapidity values 
of 0.20, 0.25, 0.35, 0.45, 0.55 and 0.65 [13,18,19]. 

From the measurements in Ref. [13], the ratio of 
data to MC simulation predictions is used to supple- 
ment the E/p measurements in Section 3 with a larger 
energy range. However, since these measurements are 
not made with the same detector, additional system- 
atic uncertainties from the test beam have to be taken 
into account [13]: 

— a fully correlated (in p and r]) 0.7% total energy scale 
uncertainty for the energies in the LAr calorimeters; 

— a fully correlated 0.5% total energy scale uncer- 
tainty for the energies in the TileCal; 

— an additional 0.4% uncertainty on the LAr unifor- 
mity for all measured energies at the same ry points; 

— an additional 1.5% uncertainty on the TilcCal uni- 
formity for all measured energies at the same r/ points; 
and 

— a 1% uncertainty from the material in front of the 
calorimeter. 
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Fig. 15: Ratio of the (E/p) measurement obtained with topological clusters to that obtained using all calorimeter 
cells in (a) the central {\r]\ < 0.6) and (b) forward (0.6 < \r]\ < 1.1) regions. The inset shows the ratio of the MC 
prediction to data. 



For single particle momenta above 400 GeV no di- 
rect measurements in a test beam or in situ exist. There- 
fore an additional uncertainty of 10% is added in quadra- 
ture to that of the 350 GeV measurement uncertainty 
in order to cover possible effects from calorimeter non- 
linearities at high energy densities and longitudinal leak- 
age [20]. 

5.2.2 Absolute calorimeter energy scale 

The absolute electromagnetic energy scale in ATLAS 
has been established using Z — ?> ee decays for the elec- 
tromagnetic LAr calorimeters and using the energy loss 
of minimum ionising muons in the TileCal. For the bulk 
of the electromagnetic LAr barrel calorimeter, the un- 
certainty on the cell energy measurement is 1.5%, and 
for the LAr presampler the uncertainty is 5% [21]. For 
the TileCal, the scale uncertainty is 3% [22] . This uncer- 
tainty does not affect charged particles with E/p mea- 
sured in situ, but needs to be considered for all other 
particles contributing to jets. 

5.2.3 Baryons and neutral hadronic particles 

Test beam measurements of protons [13,18,19,20] and 
E/p measurements for identified pions and protons (Sec- 
tion 4) have shown that the agreement between data 
and MC simulation for protons is similar to the data to 
MC agreement for charged pions. The most significant 
data to MC difference in the energy response is 25% 



for low energy anti-protons. However, on average, these 
low energy anti-protons contribute no more than 0.5% 
to the jet energy. Hence no additional uncertainty for 
charged baryons is assumed. 

In addition, the charged particle composition of the 
tracks (mainly tt , p , K ) used for the inclusive E/p 
measurements in minimum bias events and the charged 
particle composition in jets is found to be sufficiently 
similar to cause no additional systematic uncertainty. 
Thus, even if some identified particle measurements 
show differences (as for low energy anti- protons), the 
inclusive E/p measures, on average, the charged parti- 
cle response that is needed for an application to jets. 
The related uncertainty is proportional to the MC sim- 
ulation modelling of the difference in the charged track 
composition between minimum bias events and jets and 
found to be negligible. 

No test beam measurements for neutral hadronic 
particles have been carried out. Moreover, the Geant4 
models have large uncertainties. On average 10-12% 
of the jet energy is carried by neutral hadrons, mostly 
Ks, Kl and neutrons. Most of the Ks decay to pions 
before they reach the calorimeter. Hence the E/p and 
CTB measurements can be used for Ks. 

For neutrons and anti-neutrons Geant4 studies (see 
Fig. 16a and 16b) comparing alternative Geant4 hadronic 
physics models to the ATLAS-dcfault hadronic physics 
model QGSP_BERT [8] show that the (anti-)ncutron to 
(anti-)proton response ratio is determined at the 10% 
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level for particle momenta below 3 GeV and at the 5% 
level for higher momenta. Hence the (anti-)neutron re- 
sponse can be related to the sufficiently well simulated 
(anti-)proton response with these additional uncertain- 
ties. 

Few measurements are available for kaon interac- 
tion cross sections in materials. While the response of 
charged kaons is covered by the inclusive E /p measure- 
ments, the Kl response has to rely on Monte Carlo 
simulation predictions. Geant4 studies comparing al- 
ternative Geant4 hadronic physics models show un- 
certainties of ^10% with respect to the ATLAS-default 
hadronic physics model QGSP_BERT (see Fig. 16c). 
However, because of the limited availability of measure- 
ments, a more conservative uncertainty of 20% on the 
Kl calorimeter response is added in quadrature to the 
uncertainty for charged particles [23]. 



5.3 Jet energy scale uncertainty estimation 

This section combines the single particle uncertainties 
discussed in Sections 3, 5.1 and 5.2 into an expected 
shift and uncertainty on the calorimeter jet energy re- 
sponse with respect to the MC simulation prediction. 
Because of the CTB measurements used for high mo- 
mentum particles, the estimation is limited to the pseu- 
dorapidity region \r}\ < 0.8. 

Figure 17 shows the individual contributions to the 
uncertainty on the jet energy scale for \rj\ < 0.3: 
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The E/p response and CTB single particle measure- 
ments show a shift of less than 1% between data and 
Monte Carlo simulation, which propagates into an ex- 
pected shift of - 0.5-1% in the JES (Fig. 17a and 17c). 
Due to the high precision of the measurements, the un- 
certainty is small. The visible drop in the response due 
to CTB measurements is caused by a systematic differ- 
ence between the CTB data and Monte Carlo simula- 
tion for high momenta. 

The effects of the E/p acceptance and of the cluster- 
ing are small (Fig. 17b and Fig. 17e). The influence of 
the absolute electromagnetic energy scale uncertainty is 
small for low px jets, where most of the energy is car- 
ried by particles with E/p measured in situ. However, 
high pt jets have only a small fraction of energy in the 



particle momentum range of the E/p analysis and are 
hence fully affected by this uncertainty (Fig. 17d). 

The uncertainty on the calorimeter response to neu- 
tral hadrons was estimated in Section 5.2.3, resulting 
in a ~ 1% contribution to the total JES uncertainty 
(Fig. 17f). 

These individual contributions to the uncertainty 
are summarised in Fig. 18 together with the total ex- 
pected shift and uncertainty. No single component is 
dominant and depending on the jet momentum several 
components contribute at approximately the same level 
to the total uncertainty. 

Finally, the total calorimeter uncertainty on the jet 
energy scale is shown in Fig. 19 for anti-fct jets with 
R = 0.4 and R = 0.6 in the pseudorapidity range \r]\ < 
0.8. For both jet sizes the maximum expected shift in 
the jet energy scale is ^ 1% with an uncertainty of 1- 
3%. The envelope of the shift and uncertainty on the 
calorimeter JES is taken as the contribution to the total 
JES uncertainty discussed in Ref. [1]. 



5.4 Jet energy scale uncertainty correlations 

The use of pseudo-experiments for the determination 
of the JES uncertainty allows a direct extraction of the 
correlation of uncertainties between different jet mo- 
menta, pseudorapidities or algorithms by correlating 
fluctuations of different quantities within each pseudo- 
experiment. Figure 20 shows the correlation of the JES 
uncertainty between different ry and p bins. As expected, 
the correlation between neighboring bins in |?7| and px 
is almost 100%, while widely separated bins show only a 
^ 30% correlation. This remaining ~ 30% correlation is 
mostly caused by the calorimeter energy scale and neu- 
tral hadron uncertainty, which contribute identically to 
all jets. 

One specific use for the correlation of the uncer- 
tainty between jets is to compare the relative calorime- 
ter response between jets selected from two different 
categories (different reconstruction algorithm, different 
originating parton, etc.). For each jet category, the pro- 
cedure of Section 5.3 is repeated and the expected shift 
and uncertainty on the JES with respect to the MC 
simulation is derived^. For the determination of the un- 
certainty on the relative calorimeter response, the ratio 
of two jet categories is taken, cancelling out all common 
shifts and uncertainties. 

® Because the shifts and uncertainties are evaluated with 
respect to the MC simulation, any response shift already vis- 
ible at the MC level will not propagate into either the JES 
uncertainty or into the relative calorimeter response between 
jets from different categories. 
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Fig. 16: Ratio of the average calorimeter response for alternative Geant4 hadronic physics models FTF_BIC, 
FTFP.BERT and QGSP_BERT.CHIPS [8] to the ATLAS default hadronic physics model QGSP_BERT as function 
of the particle kinetic energy Skin in the range between 100 MeV and 50 GeV for (a) the ratio of neutrons n to 
protons p, (b) the ratio of anti- neutrons n to anti-protons p and (c) neutral kaons K^ ■ 



In Fig. 21a the relative calorimeter uncertainty on 
the JES between anti-fcf jets reconstructed with R = 0.4 
and R = 0.6 is shown. 

Figure 21b shows the relative calorimeter uncer- 
tainty for 6-tagged jets in tt events with respect to jets 
in inclusive di-jet events. Identified 6-jets are associated 
with a displaced secondary vertex reconstructed by the 
SVO algorithm [24], with a weight greater than 5.72. 
This weight gives an average 6-tag efficiency of 50% for 
6-jets from tt events. 

Finally, Fig. 21c shows the relative uncertainty for 
jets initiated mostly by quarks with respect to jets in 
inclusive di-jet events. Jets initiated mostly by quarks 
are selected in 7-l-jet event simulations using the se- 



lection criteria applied in the method of the direct px 
balance between 7 and jets [1]. 

Because of the small differences in particle compo- 
sition between these different jet samples, the relative 
uncertainty on the pure calorimeter response is found 
to be below 0.5% for low px and negligible for high px- 



6 Conclusions 

The average calorimeter response to isolated hadrons 
with respect to the track momentum (E/p) has been 
measured in minimum bias events at ^/s ~ 900 GeV 
and -y/s = 7 TeV. A background from neutral hadrons 
of 4-8% is subtracted from data and Monte Carlo with 
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Fig. 17: Expected shift (black markers) and uncertainty (error bars) on the relative calorimeter jet response with 
respect to the MC simulation for jets reconstructed with the anti-fct jet algorithm {R = 0.6) in the range \t]\ < 0.3 
as function of the jet transverse momentum from the uncertainty on (a) the E/p response, (b) the uncertainty 
on the E/p acceptance, (c) the uncertainty from the CTB response, (d) the uncertainty from the absolute energy 
scale, (e) the uncertainty from the clustering effect, and (f) the uncertainty from neutral hadrons. The x-axis is the 
jet transverse momentum calibrated from the EM scale to the hadronic scale using an MC-based JES calibration 
factor [1]. 
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Fig. 18: Expected total shift (black dots) and uncertainty (light blue band) on the relative calorimeter jet response 
with respect to the MC simulation for jets reconstructed with the anti-Zcf jet algorithm {R ~ 0.6) in the range 
\r]\ < 0.3 as function of the jet transverse momentum. The uncertainty components from Fig. 17 are shown 
superimposed on the expected total shift. The x-axis is the jet transverse momentum calibrated from the EM scale 
to the hadronic scale using an MC-based JES calibration factor [1]. 
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Fig. 19: Expected shift and total uncertainty on tlie relative calorimeter jet response with respect to the MC 
simulation for jets reconstructed with the anti-Zct jet algorithm (R = 0.4 and R = 0.6) in the range \t]\ < 0.3 and 
0.3 < I?;! < 0.8 as function of the jet transverse momentum. The x-axis is the jet transverse momentum calibrated 
from the EM scale to the hadronic scale using an MC-based JES calibration factor [1]. 



a systematic uncertainty of below 1%. After the back- 
ground has been removed, the agreement between data 
and Monte Carlo simulation is within 2% for particles 
with momenta up to 10 GeV and is around 5% for mo- 
menta in the 10-30 GeV range, where the statistical 
uncertainty dominates the total uncertainty. 

The calorimeter response of identified single charged 
hadrons has been measured using short-lived particle 
decays. A good agreement between data measured at 
y/s ~ 7 TcV and the Monte Carlo simulation is found 
for charged pions and protons. However, a disagree- 
ment of up to 10% is found between the Monte Carlo 
simulation and data for the difference of responses to 
low momentum anti-pions and anti-protons {{E /p)'" — 



{E/p)P). This difference is attributed to the poor mod- 
elling of the anti-proton response in the Monte Carlo 
simulation. 

The ATLAS calorimeter jet energy scale uncertainty 
has been determined for the well understood central de- 
tector region by propagating the energy response uncer- 
tainty of all particles contributing to a jet. For charged 
hadron momenta below 20 GeV, the single charged hadron 
response has been used, while for higher momenta, the 
response measured in the ATLAS combined test beam 
has been included. For the response to neutral pions 
(tt'^ — > 77), the uncertainty on the electromagnetic calorime- 
ter energy scale is dominant, while for all other neutral 
hadrons an additional uncertainty due to the limited 
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Fig. 20: Correlation coefBcient of tlic total uncertainty on the calorimeter jet response for jets reconstructed with the 
anti-kt jet algorithm (R = 0.6) in the range \ri\ < 0.8 and jet transverse momenta between 15 GeV< pt <2.5 TeV. 
The bins on the x-axis are identical to the bins on the y-axis. A dark area indicates highly correlated uncertainties, 
a light area almost uncorrelated uncertainties. 



knowledge of the calorimeter response has been incor- 
porated. 



An uncertainty of 1-3% on the response to jets in 
the calorimeter is determined for jets in the central re- 
gion of the detector with \ri\ < 0.8 and a transverse mo- 
mentum between 15 GeV and 2.5 TeV. An analysis of 
the correlation in the jet energy scale uncertainty yields 
an almost complete correlation for jets close in \r]\ oi pt, 
and a correlation of ^ 30% for jets with a large trans- 
verse momentum separation. The relative jet energy 
scale uncertainty between anti-fc^ jets of size R — 0.4 
and R = 0.6, between ^-tagged jets and inclusive jets 
in di-jet events and between mostly quark-initiated jets 
and inclusive jets in the same event sample was found 
to be small. 



A Validation of the background subtraction 

The background subtraction procedure relies on the assump- 
tion that the energy density deposited by other particles in 
the cone around a late showering track is constant as a func- 
tion of the distance with respect to the track impact point 
on the calorimeter. Although this assumption is valid for low 
track momenta, the validity of the assumption for high track 
momenta is questionable, as jet-like structures are expected 
to emerge. 

In order to address this issue, a different background sub- 
traction procedure has been developed. The background is 
still measured in events where the energy in the EM calorime- 
ter associated to the track is compatible with that of a late 
showering hadron (MlP-tagged track, in the following). In 
this case, however, the energy density of the background in 
small annuli between a cone of radius 0.1 and 0.2 is used to 
linearly extrapolate the background inside the cone of size 
0.1. The differential energy density pi = d^ E / dridcj) for an 
annulus of inner radius Rt = i x 0.025 and outer radius 
Ri+\ = (i -h 1) X 0.025 is defined as 

P« = ^ ^^ (6) 
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Fig. 21: Expected shift (black dots) and uncertainty (error bars) on the relative calorimeter jet response as function 
of the jet transverse momentum for jets in the range 0.3 < \ri\ < 0.8: (a) between jets reconstructed with the anti- 
kt jet algorithm with R = 0.4 and jets reconstructed with R — 0.6, (b) between 6-tagged jets in tt events and jets 
in inclusive di-jet events, and (c) between mostly quark-initiated jets in 7-l-jet events and jets in inclusive di-jet 
events. The x-axis is the jet transverse momentum calibrated from the EM scale to the hadronic scale using an 
MC-based JES calibration factor [1]. 



where Ai = 7r_R| . 

The dependence of pi on the distance from the track im- 
pact point on the calorimeter has been studied by making 
use of the true energy deposits as predicted by the Geant4 
simulation. The true energy deposited by the track and by 
the neutral background for late showering tracks, and by the 
neutral background in all selected events is shown in Fig. 22 
for two different track momentum bins. The energy density 
of the MlP-tagged tracks is very narrow in both momentum 
bins, as expected. The background energy density is constant 
and independent of the distance to the track impact point at 
low track momenta. This is consistent with the assumption 
of constant energy density made to evaluate the background 
with the baseline method. As the track momentum increases, 
the energy density close to the track impact point is higher. 



perhaps indicating a jet-like structure around the track. In 
both track momentum bins, the background energy density 
in the annulus 0.1 < i? < 0.2 is the same for MlP-tagged and 
all tracks. 



The same features can be observed in Fig. 23. This time, 
the reconstructed energy density is shown for data and MC 
simulation (the points corresponding to the true energy den- 
sity associated to the MlP-tagged track are kept for refer- 
ence). A comparison of Fig. 22 and Fig. 23 reveals that the 
main features observed for the true energy densities still hold 
for reconstructed energy densities: the density in the halo of 
the cone is a good measurement of the true background den- 
sity, while the core of the cone mainly contains energy from 
the MIP track. 
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The alternative estimate of the background is obtained 
by fitting the measured energy density with a linear function 
between Ri = 0.1 and Ri = 0.25 and integrating the obtained 
function between Ri = and Ri = 0.25: 



f{R,) = a + bR, 



(7) 



The obtained background estimation gives a total of 10% 
more background at high track momenta with respect to the 
baseline background subtraction procedure, consistently for 
data and MC simulation. The reason for the small difference 
is that the two background estimates differ only in the core 
of the cone, which is a region of small area. 

Since the background itself is a 10% correction to the 
(E/p) measurement, the total effect introduced by the linear 
parameterisation of the energy density is of the order of 1% 
on (E/p). Since data and MC points are affected in a similar 
way, the effect on the agreement between MC and data is 
negligible. 
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